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CN . Abstract 

I We consider the massless tricritical Ising model A^(4, 5) perturbed by the thermal oper- 

■ ator (^1^3 in a cylindrical geometry and apply integrable boundary conditions, labelled 

^r~i . by the Kac labels (r, s), that are natural off-critical perturbations of known confor- 

mal boundary conditions. We derive massless thermodynamic Bethe ansatz (TBA) 
equations for all excitations by solving, in the continuum scaling limit, the TBA func- 
^ . tional equation satisfied by the double-row transfer matrices of the lattice model 

\ of Andrews, Baxter and Forrester (ABF) in Regime IV. The resulting TBA equations 

describe the massless renormalization group flow from the tricritical to critical Ising 
^SJ . model. As in the massive case of Part I, the excitations are completely classified in 

\ terms of (m, n) systems but the string content changes by one of three mechanisms 

along the flow. Using generalized g-Vandemonde identities, we show that this leads 
to a flow from tricritical to critical Ising characters. The excited TBA equations are 
solved numerically to follow the continuous flows from the UV to the IR conformal 
Q^' fixed points. 



o 
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1 Introduction 

In integrable Quantum Field Theory (QFT), the Thermodynamic Bethe Ansatz (TBA) PIE] 
continues to be an important method for the study of massive and massless Renormalization 
Group (RG) flows including the study of excited states Moreover the Tricritical 

Ising Model (TIM), which is the simplest member A^(4, 5) of the unitary minimal series [H] 
beyond the critical Ising model 7V((3,4), remains a rich example for studying both thermal 
and boundary flows El El CD] • 

In Part I of this series (hereafter referred to as PCAI [11 J) we considered the massive 
tricritical Ising model 7V1(4, 5) perturbed by the thermal operator y^i s in a cylindrical ge- 
ometry and systematically derived the TBA equations for all excitations by using a lattice 
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approach. More specifically this was achieved by solving, in the continuum scaling limit, 
the TBA functional equation satisfied by the double-row transfer matrices of the A4 lattice 
model of Andrews, Baxter and Forrester (ABF) in Regime III J2] • In this paper we turn our 
attention to the massless tricritical Ising model which is obtained as the continuum scaling 
limit of the A4 lattice model in Regime IV. Our goal is to systematically derive and study 
the massless TBA equations which describe the renormalization group (RG) flow from the 
tricritical A^(4, 5) to critical A^(3, 4) Ising model. We apply the methods developed in PCAI 
and use the concepts and notations introduced in that paper without further elaboration. 

The layout of the paper is as follows. In Section 2 we discuss the classification of 
excitations using (m, n) systems including a description of the three mechanisms by which 
the string contents change along the flow. In Section 3 we derive in detail the massless TBA 
equations in the (r, s) = (1, 1) vacuum sector. We do not discuss the very similar derivation 
of the TBA equations in the other 5 sectors. The numerical solution of the TBA equations to 
yield continuous flows is discussed in Section 4. We finish in Section 5 with a brief discussion. 



2 Classification of Excitations and Flows 

For small perturbations, the scaling limit of the excitations in Regime IV of the A4 lattice 
model are classified by precisely the same (m, n) systems fH^ IT^ as at the conformal critical 
point and throughout the massive Regime III . Unlike the massive case, however, we 
find that in the massless regime the string content can change by one of three mechanisms 
along the flow. This was first observed empirically by direct numerical diagonalization of a 
sequence of finite-size transfer matrices approaching the scaling limit 

rnR^ liin NW , u = hlA (2.1) 

for modest values of the system size and < mR < 5 and is confirmed by our numerical 
solutions of the TBA equations. Here fi measures the perturbation strength and t is the 
departure-from-criticality variable. The mass m and continuum length scale R usually occur 
in the single combination mR. Notice that in the finite-size scaling we use the Regime III 
correlation length exponent z/ = 5/4 even though the actual jTHj correlation length exponent 
in Regime IV is z/' = 5/2. 

Let us consider the vacuum sector with boundary condition (r, s) = (1,1). The exci- 
tation energies are given by the scaling limit of the eigenvalues of the double-row transfer 
matrix D{u), or equivalently the normalized transfer matrix t{u), where u is the spectral 
parameter. The two relevant analyticity strips in the complex w-plane are respectively 

A 3A 

— < Re{u) < — , 2A < Re(n) < 4A (2.2) 

where A = 7r/5 is the crossing parameter. The excitations are classified by their string 
contents 

rrii = {number of 1-strings in strip i = 1,2} (2 3) 

Hi = {number of 2-strings in strip i = 1,2} 
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At the conformal critical point {mR = 0) the string contents satisfy the (m, n) system 

m + n = ^{Nei+Im) (2.4) 

where mi, 7712 and are even, m = {1711,1712), n = (111,112), ei = (1,0), and X is the A2 
incidence matrix with entries Ij^k = <5|j-fc|,i- For the leading excitations 1711,1712,112 are finite 
but Til ~ N/2 as N ^ 00. 

As explained in PCAI, an excitation with string content (m, n) is uniquely labelled by 
a set of quantum numbers 

/ = (/« |/(^)) = (/f), /«,..., I If ), If , . . . , /£)) (2.5) 

where the integers jj*'' G {0,1,2,...} with i = 1,2 give the number of 2-strings whose 
imaginary parts w^!^ are greater than that of the given 1-string v^^\ The 1-strings v^ ^ and 
2-strings labelled by 1 are closest to the real axis. The quantum numbers jj*'' satisfy 

>•■■>/£ >0, z = l,2. (2.6) 

For given string content [m, n) , the lowest excitation occurs when all of the 1-strings are 
further out from the real axis than all of the 2-strings. In this case all of the quantum numbers 
vanish jj*'' = 0. Bringing the location of a 1-string closer to the real axis by interchanging 
the location of the 1-string with a 2-string increments its quantum number by one unit and 
increases the energy. 

Although we do not make use of it here, we mention in passing that, for the tricritical 
Ising model, there exists a bijection between the patterns of 1- and 2-strings classifying 
the eigenvalues of the double-row transfer eigenvalues and RSOS paths on A^. This has as 
a consequence that the finitized partition functions (Virasoro characters) satisfy the same 
recursions ^2] as the one-dimensional configurational sums of the Corner Transfer Matrices 
(CTMs). 



2.1 Three mechanisms for changing string content 

Let t{u) be the normalized double row transfer matrix as in PCAI. In Regime IV the eigen- 
values t{u) of t{u) are doubly periodic meromorphic functions in the period rectangle 

period rectangle ~ ( ^ ^ {—Tcis, T^is) (2.7) 

where e is related to the departure-from-criticality variable t by 

t = -exp(-27r£). (2.8) 
However, in Regime IV the ^4 lattice model also admits the symmetry 

t{u±Tx/2 + me) =t{u) (2.9) 
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so we can restrict ourselves further to the rectangle {—^,2X) x {—7iie,7iie). In particular, 
this means that for mR > the distinction between the two strips effectively disappears 
— the two are joined along the scaling edge at Im(M) = 7re/2 in the scaling limit !^2.1\i . 
This allows zeros to move between the putative strips 1 and 2 by crossing the scaling edge. 
The combination of strips 1 and 2 in the upper half plane into one extended strip is shown 
in Figure 1. There is a complex conjugation symmetry between the upper and lower half 
planes. Note that this is consistent with the picture at the conformal point since in the limit 
mR —>■ the zeros are infinitely far below (strip 1) or infinitely far above (strip 2) the scaling 
edge as £ oo. 
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Figure 1: Upper part of the period rectangle (— ■|,2A) x {—me, me) in the complex -u-plane 
for the lattice model in Regime IV showing the putative strips 1 and 2. The schematic 
location of zeros reflects the complex conjugation symmetry and the symmetry ()2.9|1 . In 
the scaling limit the imaginary period e —>■ 00 and it is the behaviour of the zeros near the 
indicated scaling edge at Im(M) = 7ie/2 that is relevant. 

Let us consider the extended analyticity strip — | < Re('u) < Empirically, we find 
that as mR is increased from the strip 1 zeros approach the scaling edge at lm{u) = 7re/2 
from below while the strip 2 zeros approach the scaling edge from above. Potentially, this 
allows zeros to collide or migrate from one strip to the other. In fact we find that the zero 
patterns change by one of the following three mechanisms which occur at the scaling edge: 
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A. 



If 1-strings are closest to the scaling edge in both strips 1 and 2 then these 1-strings 
collide and move symmetrically in opposite directions parallel to the scaling edge until 
they reach Re(M) = 0, A. At this point the two zeros no longer contribute and they can 



be removed from the analyticity strip. This mechanism applies if I, 



(1) 



r(2) 



0. 



B. 



If a 2-string is closest to the scaling edge in strip 2 and a 1-string is closest to the scaling 
edge in strip 1 then the 2-string moves to the scaling edge and leaves the analyticity 



strip. This mechanism applies if /; 



(1) 

mi 



and Iml 



> 0. 



C. Otherwise, if a 2-string is closest to the scaling edge in strip 1, then this 2-string moves 
to the scaling edge and leaves the analyticity strip. This mechanism applies if Iml > 0. 

Similar mechanisms are also observed [Hj in the boundary flows of the tricritical Ising model. 



If 1712 = only Mechanisms B and C apply according to whether Iml = or Imi > 0. If 
mi = 1712 = then Mechanism C applies. In each of the three mechanisms exactly two 
zeros are removed from the extended analyticity strip. These three mechanisms are shown 
schematically in Figure 2. They are confirmed by the numerics discussed in Section 4. 



r(l) 








Trie/ 2 
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-AO A 2A -A A 2A -A A 2A 

ABC 
Figure 2: Schematic representation of the three mechanisms A, B, C by which the string 
contents change under the flow. In each case the two circled zeros leave the analyticity 
strip as indicated by the arrows. In mechanism A two 1-strings leave the strip whereas in 
mechanisms B and C it is a 2-string that leaves the strip. Note that only the location of the 
1-strings enter the TBA equations. 



2.2 Operator flow 

Remarkably, as we explain in this subsection, the empirical rules giving the three mechanisms 
A,B,C for changing string contents under the flow suffice to determine a map in each sector 
(r, s) between finitized characters for the UV [ttiR = 0, c = 7/10) and IR [mR = oo, 
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c = 1/2) fixed points. Since two zeros leave the analyticity strip under the flow, this map 
from the UV to IR takes the form 



l<r'<2, l<s',r<3, 1 < s < 4. 



(2.10) 



We flnd that the primary operators (r, s) of the tricritical Ising model flow to the primary 
operators (r', s') of the critical Ising model as shown in Table 1. This pattern of flows is 
consistent with the flow of operators observed with periodic boundaries but, in the presence 
of a boundary, the flow is always to the primary operator not to an associated descendant 
as can occur in the periodic case. 
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(r', s') 
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J_ 
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Table 1: Flow of primary operators (r, s) ^ (r', s') displayed in the Kac table. 



2.3 RG flow from Xi,i{q) to Xi,i{q) 

In the vacuum sector (r, s) = (1, 1) the (m, 72)^^+2 system of the tricritical Ising model 
is 

mi +ni = \{N + 2 + 1112), ■m2 + n2 = \mi (2.11) 

and these relations determine rti and n2 in terms of m\ and 1712. Similarly, let m and n be 
the number of 1- and 2-strings of the critical Ising model in the vacuum (r, s) = (1,1) sector 
satisfying the ^3 (m, n) ^ system 

m + n = \{N + m) or n = \{N-m). (2.12) 

Then m and n are given by the total number of 1-strings and 2-strings in the extended strip 
in the IR limit 

Jmi+m2-2, A [^1 + ^2, A 

m = < n = < (2.13) 

\jni+m2, B,C [^1 + 722-1, B,C. 

The mechanisms A, B, C determine which energy level flows to which energy level 
under the RG flow. As explained in ^H] and PCAI, the energies at a conformal point are 
determined by the string content m and the patterns of zeros in the complex w-plane. In 
terms of quantum numbers the precise mapping of energy levels under the flow is given by 

/ = (jf \ . . . , j« I jf ), , . . . , /(?)) ^r = (/(, . . . , rj (2.14) 
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Table 2: Mapping of the first 22 energy levels in the (r, s) = (1,1) sector under Mechanisms 
A, B or C. The energies, string contents and quantum numbers are shown in both the UV 
and IR regimes. The degeneracies of all the levels reconstruct the mapping between the 
tricritical and critical Ising characters xi,i{Q) ^ Xi,i{q)- 
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where 



A: 



B,C: 



r 

mi" 



(1) 



^2 -'ma-fc' 



mi- 



-1+k 



3 
k 



1,2,. 
1,2, 



m2+l-fc' 



A; 



1,2,. 
1,2, 



,mi - 1 
. , m2 - 1 

. . . ,mi 

,^2 



(2.15) 
(2.16) 



Details of this mapping for the first 15 energy levels in the (r, s) = (1, 1) sector are shown in 
Table 2. For given string content m, the base energy level Em is determined by the Cartan 
matrix C 

As. 



q4 



(2.17) 



The base energy occurs when the location of all of the 1-strings are further from the real 
axis than the locations of all of the 2-strings. Additional excitation energy is generated 
by permuting the order of 1-strings and 2-strings in each strip as dictated by the sum of 
quantum numbers Ij in each strip and is given by the Gaussian polynomials or g-binomials 



n I\ 



m+n 



m J g 



Ji- 



/l =0/2=0 



{(l)m{q)n' 

0, 



m, n > 
otherwise 



(2.18) 



with the g-factorials {q)m = (1 ~ Q') ■ ■ ■ (1 ~ Q'™) for m > 1 and (g)o = 1. The energy E is 
increased by one unit each time a 1-string is brought closer to the real axis by interchanging 
its location with the location of an adjacent 2-string. The product of two g-binomials is the 
generating function for the conformal spectra with given string content in each strip. The 
g-binomials satisfy the properties 



n 




' n-1 ' 


n—m 




m—l 



m rn— 1" 



fm+n] mn\m+n~\ 

V m \q - ^ L m Ji- 



(2.19) 



The character Xi i {q) is the generating function for the tricritical Ising conformal spectra 
in the (1, 1) sector. Explicitly, using the recursion ()2.19|) . we decompose it into three terms 



4,Af+2/ N 1- i 



mCm 



m\+n\ 




m.2+n2 


mi 




rn2 



2I0 j2 



g|(™l~™l"^2+mi) 



mi+ni— 1 
mi — 1 



mi ,m2 even 

I „m2 fmi+ni-l 
"I" y [ mi-1 



m2+n2— 1 
m2 



+ q 



mi 



■m2+n2-l 
1712 — 1 



mi+rti — 1 
mi 



(2.20) 



m2+n2 
m2 



These three terms correspond precisely to the energy levels effected by mechanisms A, B 
and C respectively. So simply reading off the conformal energies from the respective zero 
patterns after applying each mechanism we find the following mapping between finitized 
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characters 



m i ,m2 even 

j(n2 — l)mi + (n2 — l)m2 fmi+ni — 1 
[ mi-1 



J „n2(mi — l)+n2(r7i2 — 1) 



mi+ni— 1 




m2+n2— 1 


mi — 1 




m2 — 1 



"12 +"2-1 
m2 



.n2mi+n2m2 



mi+jii — 1 
mi 



g-i J] ^ |gn2(mi-l 



jni+rti — 1 
mi — 1 



m2+"2— 1 
m2 — 1 



m even mi even 



+ q 



("2- 



-l)mi mi+ni — 1 
[ mi-1 



m2+n2-l 
m2 



+ q 



n2m\ 



mi+ni— 1 
mi 



m2+n2 
m2 



m2+n2 
m2 



(2.21) 



q q^^"^^ q^"^^ """" \ q" 



m even mi even 



m^— mmi J m— ^mi+2 



N + m — mi +2 
2 

mi — 1 



mi —2 
2 

m— mi+1 



g 48 \ ^2*"^ 



JV + m. 
2 

m 



+ q-'^' 

3,N 



Af + m — mi 
2 

mi — 1 



+ 



iV+m — mi 
2 

mi 



2 

m— mi 



Xia [q) 



m even 



Notice that, after the mapping, the g-binomials of strip 2 are with respect to 1/q. This is 
because strip 2 is turned upside down when it is placed on top of strip 1 to form the extended 
strip. These g-binomials are naturally replaced by g-binomials in q using ()2.19|1 . All integers 
are then eliminated in favour of m and mi using the appropriate relations which apply to 
the mechanism corresponding to each of the three terms. The final equality holds because 
of the remarkable generalized g-Vandermonde identity which is proved in the Appendix 



N+m 
2 

m 



E 



^|mj— mmi J ^m— |mi+2 



JV + m — mj +2 
2 

mi — 1 



mi —2 



m— mi +1 



(2.22) 



+ q 



-mi 



AT+m — mi 
2 

mi — 1 



m—nii 



+ 



A^ + m — mi 
2 

mi 



2 

m— mi 



The finitized Ising character Xi'iiq) is not the usual finitized character. In the limit 
g — > 1 it gives the correct counting 



i,mx?:fw=Em = ['in„ 



(2.23) 



m even 



whereas the usual finitized Ising character gives A3 counting 



limg « q 



im^ 



E(S) = [^3"]m 



(2.24) 



where and A3 denote the adjacency matrices. 
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2.4 RG flow in other sectors 



The analysis of the flow using the three mechanisms A, B and C can be extended to each of the 
sectors (r, s) with r = l,2,s = l,2,3. In this way we obtain six generalized g-Vandermonde 
identities as follows: 



Xi,i N even, m even: 



(7V+m)/2 
m 



E 



mi even 



3mf /2— mmi 



{N+m-mi)/2 
mi 



mi/2 
m— mi 



+ ^-mi 
X3,l 



(Ar+m-mi)/2 
mi — 1 



r(mi-2)/2 
I m— mi 



+ q 



m-7mi/2+2 



(Ar+m-mi+2)/2 
mi— 1 



(mi-2)/2 
m— mi+1 



even, m odd: 

mi even 



(Af+m+l)/2 
m 



3m^/2— mmi 



{[ 



(iV+m-mi+l)/2 
mi 



mi/2 
m— mi 



+ q 



-mi 



(Af+m-mi+l)/2l r(mi-2)/2 



mi — 1 



m— mi 



+ q 



.m-7mi/2+2 



(Af+m-mi+3)/2 
mi— 1 



(mi-2)/2 
m— mi+1 



(2.25) 



(2.26) 



X2,2 N even, m even: 



+ q 



{N+m)/2 
m 



-m+5mi/2+l 



^3mf/2-mmi | 

mi even 



(W+m-mi)/2 
mi 



mi/2 
m— mi 



(Ar+m-mi-2)/2 
mi 



mi/2 
m— mi— 1 



+ q 



m— 5mi/2+l 



(Ar+m-mi)/2 
mi— 1 



mi/2 
m— mi+l 



(2.27) 



X2,i N odd, m even: 



(Af-|-m-l)/2 
m 



^3mf/2-mmi |^mi/2 

jTti odd 



(Af+m-mi-2)/2 
mi 



(mi+l)/2 
m— mi 



+ q 



-mi/2 



(Af+m-mi-2)/2 
mi — 1 



Rmi-l 



)/2 
mi 



,m— 3mi-|-3/2 



(Ar-|-m-mi)/2 
mi — 1 



(m 
m 



i-l)/2ll 
-mi+lj J 



(2.28) 



Xi,2 A" odd, m even: 



(A?"+m-l)/2 



^ ^ ^3mJ/2— mmi ^^m— 3mi-|-3/2 



(Ar+m-mi)/2 
mi — 1 



mi odd 



+ q 



-m+2mi+l/2 



(Ar+m-mi-2)/2 
mi 



(mi-l)/2l ^-mi/2 
m— mi — 1 J 



q~ 



X3,2 A^ odd, m odd: 



+ q 



(N+m)/2 
m 



-m+2mi+l/2 



3m|/2— mmi J m— 3mi-|-3/2 



{N+m-mi)/2 
mi 



(iV-|-m-mi-|-l)/2 
mi— 1 



(mi-l)/2l 
m— mi+l J 

(mi-l)/2l 1 
m— mi J J 



(mi-l)/2 
m— mi+1 



mi odd 

(JV+m-mi-l)/2 
mi 



(mi-l)/2l , -mi/2 
m-mi-lj ~ y 



(Ar+m-mi-|-l)/2 
mi 



(mi-l)/2 
m— mi 



(2.29) 



(2.30) 



These identities are simplified and proved in the Appendix. 
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3 TBA Equations in Regime IV 

Recall from PCAI that the normalized double row transfer matrix is defined by 



where 

and 
with 



t{u) = Sr,s{u)S{u) 



S{u) 



^9i(u + 3/,p)?9i(m + Z,p)J 



Diu) 



Sr.Ju) 



hriu) 

hs{u) 



M'^u-3l,p)^i{2u + l,p) 

-iyhr{u - ^L)h-r{u + ^L)hs{u - ^R)h_s{u + ^r) 

^i{l,p)^i{u + (3 - r)l,p)^i{u + (1 - r)l,p) 

'di{u,p)'di{u — l,p)'di{u + 21, p) 
1^4(1, p)Mu + (3 - s)l,p)Mu + (1 - s)l,p) 



(3.1) 

(3.2) 
(3.3) 

(3.4) 
(3.5) 



di{u,p)'di{u — l,p)di{u + 21, p) 
Moreover, the normalized transfer matrix satisfies ^Hj the universal TBA functional equation 

t{u)t{u + I) = I + t{u + (3.6) 

independent of the boundary condition (r, s). For Regime IV, the nome is pure imaginary 
with p = ie~'^^. 



3.1 UV Massless TBA: (r, s) = (1, 1) 

In this section we derive the TBA equations for the (r, s) = (1,1) boundary by solving the 
TBA functional equations in the scaling limit for even N. We follow closely the derivations 
in ^SIE] and PCAI TTJ. The derivation for other boundary conditions is similar. We begin 
by factorizing the eigenvalue t{u) of t{u) for large as 

tiu) = f{u)g{u)l{u) (3.7) 

where f{u) accounts for the bulk order- behaviour, g{u) the order-1 boundary contributions 
and l{u) is the order- 1/A^ finite-size correction. We will solve for /(«), g{u) and then l{u) 
sequentially. 

For the order- A^ behaviour the second term on the RHS of the TBA functional equation 
()3.6p can be neglected giving the inversion relation 

/(n)/(n + /) = l. (3.8) 

In the physical strip 1, the solution [20J with the required analyticity is 





TT 
4 


5m 
2 ' 


|t|5/2)n 


2N 


u 


1 + 


5m 
2 ' 


|t|5/2)J 





(3.9) 
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The solution in strip 2 satisfies the same inversion relation and is given by the symmetry 

(E3D. 

In the two analyticity strips labelled by j = 1, 2 we define generically for the functions 
h = t, f, g,l the notations 

hi{x) = h{- + —), |Im(x)|<7r (3.10a) 
^ o 

tx 

hJx) =h(3l + —), |Im(x)|<7r (3.10b) 
5 

Hi{x) =l + h{x), H2{x) = l + h2{x). (3.10c) 
and we assume the relevant functions have the scaling form 

/;,i(x) = lim /i(a; + logA^) (3.11) 

A'"— >oo 

For example, we see that 

log/i(x) = -4/xV (3.12) 

As in Regime III, we next have to solve by Fourier series the inversion relations for the 
order- 1 boundary terms 

TT TT 

9i{x-i-)gi{x + i-) = 1 (3.13a) 

TT TT 

92{x -i-)g2{x + i-) = Gi{x). (3.13b) 

To find gi{x) we need to consider the zeros and poles introduced by the prefactor in ()3.H) . 
We find the order-1 zeros of D{u) cancel exactly the poles of Si^i{u). Taking into account 
periodicity, the prefactor exhibits poles at m = — | + ame, ^ + crivre, 21 + ame, 4/ + ame, 
{±if, f ± «f , ^ ± if, ^±if and double zeros at m = f + icrvre, 3/ + iavre, f ± if, 
^ ± where a = 0, ±1. The solution for gi{x) with this analyticity but restricted to the 
strip |Im(a;)| < 37r/4 is 



"3\2 8' 




t|5/4) 




t|5/4) 


(/4l 2 8' 




t|5/4) 




t|5/4) 



(3.14) 



By comparing the expected pattern of zeros and poles of g2{x — iTT/2)g2{x + 27r/2) with 
the analyticity of Gi{x), we find that this functional relation is satisfied inside the strip 
|Im(x)| < We observe that g2{x) / gi{x) is free of zeros and poles in this strip. Similarly 
Gi{x) is analytic and nonzero, but only inside a narrower strip |Im(x)| < ^. Hence in this 
smaller strip 

log5(2(a;) = loggi{x) + e* logG'i(x) (3.15) 
where the kernel in Regime IV is given by 

= o,.o /.^ • (3-16) 
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Note that the signs of gi{x) are choosen to match the corresponding expressions in the 
tricritical hmit t — 0. They can be determined numerically from the eigenvalues of the 
transfer matrix with real u in the physical range < u < 1/2. 
The functional equations for the finite-size corrections 

h{x - j)h{x + j) = T2{x) (3.17a) 

can be converted to Nonlinear Integral Equations (NLIE) by standard techniques 122] 
where the key input is the analyticity determined by the patterns of zeros. Suppose that 1- 
strings are located in the upper half plane in the extended strip 1 at {//2 + it>i, l/2 + iVmi} 
below the scaling edge and at {1/2 + iwi, //2 + iwmi} above the scaling edge. Then the 
zeros in strip 2 are determined by the symmetry ()2.9|) and occur at {3/ + i^iie — Vi), 31 + 
i^ire — Vmi)} and {3/ + i{7ie — Wi), 31 + i(7TS — Wmi)}- To account for these zeros, we note 
that the functions 

^ ^ ^2(^f - mu^l + \t\^) fT Mi-2 - \t\'^)Mq + .3 

'^""^ y ^i(^f - ^^Jt|i)^i(^f + M ^i(^f - z^, |t|t)^i(zf + z^, |t|i) 

^ . . fr ^3(^f - |t|i)^3(^f + \t\l) - \t\'^)Mq + I'^l^) .3 

are free of zeros and poles inside their respective analyticity strips. The products of elliptic 
functions satisfy the inversion relations li{x — '^)li{x + y) = 1, i = 1,2 and are doubly- 
periodic. 

It follows that 

logti(a;) = log/i(x) + log(?i(x)+e*logT2(x) + Ci (3.20) 

^y^^^ ^i(^f-^^,|t|i)^i(^f + ^^,|t|i) ^ ^i(^f-^^,|t|i)^i(zf + ^^,|t|i) 

logt2(a;) = log/2(x) +log^2(a;) +e*logTi(s) -£*logG'i(x) + C2 (3.21) 

™i ^4(^1 - |t|5)^?4(^f + |t|i) ^ , ^4(^1 - |t|t)^4(^f + |t|t) 



+ / log 5 ^ h > lo£ 



,=1 ^3(^f-^^,|t|^)^3(^f +^^,|t|^) ' fcr"°^3(^f-^^,W^)^3(^f 

where £(x) is the kernel ()3.16|) . As in |^15j, the integration constants Ci vanish. So now 
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taking the scaling limit yields 

logti(a;) : 



+ log 



+ k*logf2{x) 



jie^ — -\/2 + ji'^e"^ 
-{(5f^ -X- log/x)j + ^ logtanh ^ -x- log/x) 



+ log tanh 

3=1 " fe=l 

log4(x) - -4e-^ + log ^"""^ ^ ^ ^ ^ V >K logri(x) 



(3.22) 



+ ^ log tanh - + x + log 

i=i 



/i) j + log tanh 

k=i 



+ x + \og^i) (3.23) 



where the kernel is 



k{x) 



2% cosh x 

In deriving this result we have assumed that the zeros scale as 



Vj ~ 



\og\t\ , Pf^ _7rs , 



4 
log 



+ 



4 5 2 5 



or, more precisely, the scaled locations of the zeros are defined by 

= lim (5^. + ^logl^l) = lim f5'u,- + log 4) 
= lim f5wfe + ^logltl) = lim f 5wfc + log 



ii 

N 

Close to the UV fixed point, when mR is small, we must have (3f^ < and (3^^^ > 0. 
If we define the rapidity •& and pseudo-energies ei{'d) by 

e-^'(^) = lim uU - log-^) = U{d - log/x), 

where // = mR/A we obtain the TBA equations 



(3.24) 

(3.25) 
(3.26) 

(3.27) 
(3.28) 



i = 1,2 



(3.29) 



ei(i?) = mRe^ - log 



mi 



1 

2^ 



-V2- 



- ,log(l + e2f2) 
cosh(i? - 1?0 



^logtanh - -) - ^logtanh ( 



(2) 



) 



e2(t?) = mRe-"^ - log 



e^-V2 + e-^ 27r 



2 2' 

,log(l + e--^(^')) 
cosh(^? - -(?') 



(3.30) 



- 5^ log tanh + ) - J] log tanh + ^) 



fe=i 
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To find the locations (3f^ , jSj^^ of tlie 1-strings consider tlie functional equations 

TT TT 
tl{x - l-)ti{x + t-) = l + t2{x) 

t2ix - t^Mx + t^) = l + h{x) (3.31) 

at X = —TTi/2 + 5vj, X = —TTi/2 + 5wk respectively. Since the right-hand sides vanish, in the 
scaling limit this implies 



t2(/?r^---log/i) = -l = e-^ -\ J = 1,2,..., mi 



,(1) TTz , , ^ _„m, 
^ 2 

,(2) vri , , ^ _.„{2), 



tm ' - y - log/i) = -1 = e-"'^ A; = 1, 2, . . . , m2 (3.32) 



or 



e2(/?f -it) = ^i^'^ J = 1,2,..., mi 



eil/^r-y) = k = l,2,...,m2 (3.33) 

where n^l\ are odd integers. These integers are given by their values ^H] in the UV 
limit as determined by winding numbers 

nf^ = 2(mi -j)-m2 + l + 2lf\ j = 1, 2, . . . , mi 

= 2(m2 - A;) - mi + 1 + 2/f \ A; = 1, 2, . . . , ms (3.34) 

where / = . . . , Imi\Ii'\ ■ ■ ■ , Iml) are the quantum numbers ()2.5|) . These integers nf'' 
can change during the flow due to winding of phases. For numerical purposes, a more useful 
form of these auxiliary equations is obtained by replacing d with Z?® — y in the TBA 
equations ()3.3()|1 . Similar equations can obtained for the locations of the 2-strings. 

Repeating the same calculation as in Regime III, we find the finite-size corrections to 
the scaled energies in Regime IV are 

1 , ^ , s RE(R) coshx , 1 , , , 

-logDi(x) = ''-^ + o(-) (3.35) 

with 

mi 1712 TD roo 

RE{R) = mR\Y}-^''' - ^ / d^e'' log(l+e-^iW) (3.36) 
i=i fc=i 

In the UV limit niR — > we recover the critical TBA equations of jT3]. Explicitly, 
setting 'd = X + log /x, we find that as /x — the location of the 1-strings scale as 

= yf^ + log /i, i = 1, 2, . . . , mi 
/3f = -yf-log/i, A; = l,2,...,m2 (3.37) 
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It follows that in the limit mR the finite-size corrections to the scaled energies are given 
exactly by 

7 1 



1 27r 

-\ogDi{u) = — sinSu 



2 rui 



240 4 



i=i j=i 



(3.38) 



with the finitized partition function 



\mCm 



mi , m2 even 



mi + rii 




m2 + n2 


mi 


q 


1712 



(3.39) 



where the modular parameter is g = exp(— Trsin SwM/iV) for M double rows. 

We have chosen to write the TBA equations in terms of both strips 1 and 2. However, 
in accord with the existence of a single extended strip, the set of TBA equations ()3.3()j) . 
()3.33p . ()3.34|) and ()3.36|) can be written in terms of a single strip by using the symmetry 
^2(1^) = ei(~''^)- Explicitly, extending strip 1, we obtain the TBA equations 



61 (^) 



mRe^ — log 



-Z]logtanh(^--) 



, log(l + e 



cosh('(9 + 



1=1 



RE{R) = mR^e-^- 



(1) mR 



2n 



/oo 
d^e^ log(l + e-^i('')) 
-00 



where m = mi + m2 and 



PI 



(2) 



k = 1, 2, . . . , m2 



The auxiliary equations become 



(1) • 
TTZ, 



1,2, 



, m 



where 



n 



(1) 



m - 2£ + 1 + 2J^ 



1,2 m 



and the quantum numbers le in the extended strip are given by 



n2 + l 
n2 



(1) 
^(1) 



1,2, . . . ,mi 
mi + 1, mi + 2, 



, m 



(3.40) 
(3.41) 

(3.42) 
(3.43) 
(3.44) 

(3.45) 



with n2 = mi/2 — m2. This symmetry between strips 1 and 2 is manifestly broken in the UV 
limit mR and the IR limit mR 00. Also these TBA equations need to be modified 
in the intermediate regime for Mechanism A levels after collision of the 1-strings. 
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3.2 IR Massless TBA: (r, s) = (1, 1) 

For large rriR we work with the extended strip 1. In this regime the total number of 1-strings 
is either m = mi + m2 for Mechanism B, C or m = mi + m2 — 2 for Mechanism A with 
quantum numbers /' = {I[,l2, . . . , I'^) given by (j2.14j] to ()2.1(j|l . Setting 

e'{^) = ei{^) + mm = e2(-'^) + mm (3.46) 

we now obtain the single TBA equation 



e'(t9) = mRe"^ - log ^ ^ ^ o x , 



1? 



+ 72 + 6-^ 1 p^^,log(l + e-W) 



e 



y/2 + 2tt cosh(79 + ^?') 



-J]logtanh(^-^) (3.47) 

i=i 

with scaling energies 

T~) /*00 

i?£;(i?) = mi? y e-^^ - ^ / rf^? e'' log (1 + e'^'^'')) (3.48) 

The auxiliary equations are 

e'(y-/5;) = n[m, £=l,2,...,m (3.49) 

where 

n[ = 2(m - ^) + 1 + 2/;, £ = l,2,...,m (3.50) 



The difference between this equation and ()3.44|1 . as reflected in fl3.46|l . arises because of the 
mni winding between the previous position of the scaling edge at Tiie/2 (reference point for 
■(9 = 0) and the new position of the scaling edge at nie. 

Setting 'd = X + log/x in the IR limit mR — > oo, the locations of the 1-strings scale as 

/?; = y; + log/i, ^=1,2,. ..,m (3.51) 

and the above equations reduce to the TBA equations of the critical Ising model. Indeed, 
the pseudo energies decouple giving the energy of the usual massless free fermions 

e'(^?) ~ mRe-^ (3.52) 

so that the auxilliary equations ()3.49|) become 

4e-^^ = -[2(m-£) + l + 2/;]7r (3.53) 

In this limit the finite size corrections for the scaled energies are 



1 , / N TT cosh X 
- log Di[x 



48 . 



2 O .V / ^ 

with the partition function 

4r (g) = q-'"' Yl 1 



E' = ^ + J2l', (3.54) 



r N+m-2 ' 
2 

m 



(3.55) 
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4 Massless Numerics 



The TBA equations of the previous section can be solved numerically by an iterative pro- 
cedure. There are, however, some subtleties. The process starts with initial guesses for the 
pseudoenergies ei{'d) and 1-string locations, close to the UV or IR fixed points. The flow 
is followed by progressively incrementing or decrementing mR. At each value of mR, the 
TBA equations are used to update the pseudoenergies ei{{}) and then these are used in the 
auxiliary equations to update the locations of the 1-strings, and so on, until a stable solution 
is reached. Typically, the UV form of the equations is stable for small values of mR, the IR 
form is stable for large values of mR and there is an intermediate range of values of mR for 
which both forms are stable and converge (with a precision of five decimals places) to the 
same values for the scaled energies and the locations of each of the 1-strings. In all cases 
these numerical flows confirm the three mechanisms A, B, C. 

One numerical difficulty is related to the determination of the location of the 1-strings 
in strip 2. This problem arises because in the UV form of the equations cannot be 
obtained by direct iteration of the auxiliary equation. This problem is solved by inverting 
the phases 

log tanh ( ^^ +-) (4.1) 

A more serious problem relates to Mechanism A cases for which two 1-strings collide to form 
a (short) 2-string with complex coordinates 

/?S'/3S^/5±^7 (4.2) 

For these cases there is an intermediate range of values of mR which requires a modification 
of the TBA equations. Such short 2-strings were studied [22] in the context of the Yang-Lee 
scaling theory. Unfortunately, due to instabilities in our equations, we have been unable to 
numerically solve the Mechanism A equations throughout the intermediate region. 

We present some typical numerical results in a series of figures. In Figure 3, we compare 
the groundstate scaling energy in the (r, s) = (1,1) sector with the scaling energy for periodic 
boundary conditions. In Figures 4 and 9 we show the fiow of scaling energies in the sectors 
(r, s) = (1, 1) and (r, s) = (3, 1) respectively. A dashed curve is used to guide the eye in the 
intermediate regime of the Mechanism A level. In Figures 5 and 6, we show the fiow of the 
scaling energy and 1-strings for the Mechanism A level in this sector with string contents 
(mi,m2) = (4,2) and UV quantum numbers / = (0, 0, 0, 0|0, 0). The dashed curves in the 
intermediate regime are schematic and have not been calculated from the solution of the 
TBA equations. For comparison, we show in Figures 7 and 8, the fiow of the scaling energy 
and 1-strings for the Mechanism B level with string contents (mi, m2) = (6, 2) and quantum 
numbers I = (0, 0, 0, 0, 0, 0|1, 1). The fiow of the scaling energies and 1-strings for arbitrary 
Mechanism B and C levels can be calculated throughout the fiow by numerical solution of 
the TBA equations. The mechanism A levels can be calculated right up to the point where 
the two 1-strings collide. Note the linear regimes in the UV and IR for the fiows of 1-strings. 
This corresponds to the assumed limiting scaling of the locations of these 1-strings. 
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Figure 3: Flow of groundstate energy Cggp = —12RE{R)/'k versus \og{mR) for periodic and 
(r, s) = (1, 1) boundary conditions. The difference arises from the boundary term in the TBA 
equations. The energy is a decreasing function of mR for periodic boundary conditions by 
Zamolodchikov's c-theorem but this theorem does not apply with fixed boundary conditions. 

5 Discussion 

In this paper we have used a lattice approach to derive TBA equations for all excitations 
in the massless renormalization group flow from the tricritical to critical Ising model. The 
excitations are classified according to string content which changes by one of three mecha- 
nisms A,B,C along the flow and leads to a mapping between finitized Virasoro characters. 
With the exception of the intermediate regime for Mechanism A flows, the TBA equations 
can be solved numerically by iteration. It would be of interest to compare our results with 
the results of the Truncated Conformal Space Approximation. 

Although, the tricritical Ising model is superconformal, the boundary conditions applied 
in this paper break the superconformal symmetry. It would be of interest to investigate the 
pattern of the superconformal flows between fixed points corresponding to superconformal 
boundary conditions [211 E] • It would also be of interest to extend the analysis of this paper 
to the complete flow for periodic boundary conditions. 
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Figure 4: Flow Xii(Q') ^ Xiiil) of scaling energies — Cgff/24 = RE{R) /2ti versus \og{mR) 
in the (r, s) = (1,1) sector. The degeneracies of the levels are shown in the margins. The 
intermediate region of the Mechanism A levels (shown dashed) are schematic and have not 
been obtained from numerical solution of the TBA equations. 
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Figure 5: The scaling energy — Cgff/24 = RE{R)/2tt versus log (mi?) for the Mechanism A 
level in the (r, s) = (1,1) sector with string contents (mi, 7712) = (4,2) and UV quantum 
numbers / = (0, 0, 0, 0|0, 0). 
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Figure 6: Flow of six 1-strings versus log(mi?) for the Mechanism A level in the (r, s) = (1, 1) 
sector with string contents (mi, 1712) = (4, 2) and UV quantum numbers / = (0, 0, 0, 0|0, 0). 
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Figure 7: The scaling energy — Cgg-/24 = RE{R)/2tt versus log(mi?) for the Mechanism B 
level in the (r, s) = (1,1) sector with string contents (mi, 7712) = (6,2) and UV quantum 
numbers I = (0, 0, 0, 0, 0, 0|1, 1). 




Figure 8: Flow of eight 1-strings versus log (mi?) for the Mechanism B level in the 
(r, s) = (1,1) sector with string contents (mi,m2) = (6,2) and UV quantum numbers 
I = (0,0,0,0,0,0|1,1). 
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Figure 9: Flow xiiil) ^ Xisil) of scaling energies — Cgff/24 = RE{R)/27i versus log(mi?) 
in the (r, s) = (3, 1) sector. The degeneracies of the levels are shown in the margins. The 
intermediate region of the Mechanism A levels (shown dashed) are schematic and have not 
been obtained from numerical solution of the TBA equations. 
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Appendix 

In this Appendix we prove the generahzed g-Vandermonde identities of Section 2.4. We sim- 
phfy the identities to show they are special cases of general identities obtained by Bender pH] . 
We then give an elementary proof of these identities using induction. 

To simplify the identities of Section 2.4 we set mi = 2k or mi = 2k + 1 depending on 
the parity of mi and set n = {N + m)/2 or n = (A^ + m ± l)/2 as appropriate. This reduces 
the six identities to four identities 
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After some recasting, a surprising mod 3 property emerges in the terms of these identities 
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Setting (. = 3k, 3k + 1,3k + 2 mod 3 reduces the four identities to just two identities 
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For (n, m) 7^ (0, 0) there is an additional identity 
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These are in fact special cases^ of identities due to Bender 
Generalized g-Vandermonde Identities: 
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Proof: For n = m we have £ = m so LHS = RHS = 1 for — m < a < 2m + 1. We now 
proceed by induction on n. Suppose that — m < a < 2?7, + 1 and — (m — 1) < a + 1 < 2n + 1, 
that is, —m < a < 2n. Then 



n+l- 



n 1 I ^n—m+1 



£ 



m— 1 



n] _|_ (?i-^)+{^-m+l) 



L(2^+2-a)/3j (e-rn) 



n-[(£+l+a)/3j 



m— 1 



L(<?+a)/3j 



n-[(£+2+a)/3j 



L(<?+l+a)/3j 



'^^l{2i+2-a)/3](i-m) 



[(2£+2-a)/3j (^-m) 



n-L(^+l+a)/3j 



n+l-[(£+l+a)/3j 



L(^+a)/3j 
l—m 



n-[(£+l+a)/3j 



L(£+a)/3j 



Now suppose that a = 2n + 1. Then only the terms i = n and i = n + l survive on the RHS 
and 



[(2£+2-a)/3j {£~m) 



n+l-L(^+l+a)/3j 

n+i-e 



L(£+a)/3j 



^^L(2^-2n+l)/3j(£-m) 



n+l-[(£+27i+2)/3j 

n+i-e 



[(£+2n+l)/3j 

e-m 



n ~\ I n—m+1 



n 
m—1 



rn+11 
L m J 
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